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Abstract 

Lunar laser ranging (LLR) has been a workhorse for testing general relativity over the past 
four decades. The three retroreflector arrays put on the Moon by the Apollo astronauts and the 
French built array on the second Soviet Lunokhod rover continue to be useful targets, and have 
provided the most stringent tests of the Strong Equivalence Principle and tire time variation of 
Newton’s gravitational constant. The relatively new ranging system at the Apache Point 3.5 
meter telescope now routinely makes millimeter level range measurements. Incredibly, it has 
taken 40 years for ground station technology to advance to the point where characteristics of 
the lunar retroreflectors are limiting the precision of the range measurements. In this article, 
we review the gravitational science and technology of lunar laser ranging and discuss prospects 
for the future. 



1 Introduction 


On .July 20. 1969. humans landed on the Moon for the first time. The following day, the Apollo 11 
astronauts deployed the first Laser Ranging Retrorefiector (LRRR or LR 3 ). Within a month, return 
photons were successfully detected at several observatories [7], Two more retrorefiector arrays were 
placed on the Moon during the Apollo 14 and 15 missions. In addition, two French arrays were on 
Russian Lunokhod rovers carried on Luna landers, however, only the array on Lunokhod 2 has been 
a viable target. These four retrorefiector arrays remain visible today and the laser ranging data 
collected over the past 40 years has dramatically and continually increased our understanding of 
gravitational physics along with Earth and Moon geophysics, geodesy, and dynamics. They remain 
the only operating lunar experiment from the Apollo era. and more incredibly, only recently has 
the ground station technology advanced to the point where uncertainties associated with the lunar 
arrays are limiting the range measurement precision. 

The Apollo 11 and 14 arrays consist of 100 fused silica circular opening cubes (3.8 cin diameter 
each) with a total estimated lidar cross-section (based on the intensity of the diffraction pattern 
of the array at the position of the receiver in the far field) of 0.5 billion square meters. Apollo 15 
has 300 of these cubes and therefore about 3 times the lidar cross-section and is the lunar array 
with the highest response. The diffraction pattern of the Apollo arrays has a bright -central iobe 
at the center with six surrounding lobes with an effective tophat function 8 arcsec in diameter for 
small incidence angles [3, 381. They produce an approximately 20 km diameter spot on Earth [231 • 
This spread is sufficient to cover the velocity aberration due to the Moon’s motion, so the cube's 
reflective face angles were not intentionally spoiled (deviation from 90 degrees). 



Figure l: From left to right, the Apollo 11, 14, and 15 retrorefiector arrays. 

The two Lunokhod arrays consist of 14 triangular shaped cubes, each side 11cm. Shortly after 
landing, the Lunokhod 1 array ceased to be a viable target - no ground stations have since been 
able to get returns from it. It is also very difficult to get returns from Lunokhod 2 during the lunar 
day. The silver rear coating and larger size of the Lunokhod cubes makes them less thermally 
stable which dramatically reduces the optical performance when sunlit. 

Lunar laser ranging (LLR) is performed by measuring the round trip light travel time between 
a ground transmitter and the retrorefiector. Early LLR measurements had a precision of about 20 
cm. Since 1969, multiple stations have successfully ranged to the lunar retroreflectors. however, 
two stations have dominated LLR data generation: McDonald Laser Ranging System (MLRS) in 
Texas (since 1969) [62] and Observatoire de la Cote d’Azur (OCA) in Grasse, France (since 1985) 
[56]. The vast majority of their lunar data comes from the array with the highest lidar cross section: 
Apollo 15. These stations have increased the range precision by a factor of 10 over the years to 



Figure 2: The Russian Lunokhod 1 lander with the French built retroreflector array sticking out 
the left side. 

the level of about two centimeters. Recently, the Apache Point Observatory Lunar Laser-ranging 
Operation (APOLLO) has begun contributing high-quality data at the millimeter level [5, 42, 39]. 

Poor detection rates are a major limiting factor in past LLR. Taking into account velocity 
aberrations, optical performance of the ground station, and other systematic effects, the overall 
round-trip loss is of order 10~ 21 , mostly due to the r 4 loss from the Earth-Moon distance. Because 
of this heavy loss of light, not every laser pulse sent, to the Moon results in a detected return 
photon, leading to poor measurement statistics. MLRS typically collects less than 100 photons 
per range measurement with a scatter of about 2 cm. However, the new APOLLO instrument at 
Apache Point has overcome this limitation. 

The large collecting area of the Apache Point 3.3 m diameter telescope and the efficient 
avalanche photodiode arrays used in APOLLO result in thousands of detections (even multiple 
detections per pulse) leading to a statistical uncertainty of about 1 mm. The dominant random 
uncertainty per photon in APOLLO stems from the orientation of the reflector array and the as- 
sociated spread of pulse return times [5]. Additionally, systematic errors associated with lunar 
arrays, such as regolith motion and thermal expansion of the array, start to become significant at 
the millimeter level of precision. 

Each ground station records the single-photon reflection events, which are then combined into 
normal points that are adjusted for station-specific corrections. A typical normal point is generated 
from 5 to 20 minutes of ranging. The normal points are then submitted to a central archive with 
the International Laser Ranging Service [55] which makes them available to the public. Auxiliary 
measurements, such as environmental conditions (temperature, pressure, etc.) are also recorded, 
as these are required to further correct the data for atmospheric effects. A detailed model of the 
solar system ephemeris is then used to perform a least-squares analysis to estimate various model 
parameters from the measured data [75]. 

Arguably, the most fruitful analysis of LLR data is for tests of General Relativity (GR). Rela- 
tivistic effects typically show up at particular frequencies (such as the synodic frequency) making 
them separable from most other parameters. These frequency signatures make the analysis possible 
given the plethora of other effects with large uncertainties. In the following sections, we summarize 
the main tests of GR performed with LLR. These include tests of the Equivalence Principle, the 
variation of the gravitational constant, the inverse square law, and a preferred-frame. We then 
discuss the current state of LLR data analysis. Finally, we describe the next generation of lunar 
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Figure 3: Improvements in the ground station technology over the past 40 years have increased I he 
range precision by 2 orders of magnitude. Uncertainties associated with the existing retroreflectors 
are now becoming a limitation. 

laser ranging instruments and discuss the possibilities that the future of lunar and planetary laser 
ranging may hold. 


2 Equivalence Principle 

The gravitational acceleration of massive bodies toward other bodies is dependent on the non-linear 
properties of gravity within metric theories of gravity [48]. Tracking this acceleration provides a 
measurement of how gravity pulls on the gravitational binding energy and how gravitational binding 
energy affects inertia. This probe of non-linear gravity is explicitly singled out in measurements of 
the Parameterized Post-Newtonian parameter 0 discussed below, but, it is also implicitly contained 
within the Einstein Equivalence Principle. 

The Equivalence Principle (EP), which states the equality of gravitational and inertial mass, 
is central to the theory of General Relativity. The EP comes in two flavors: the weak (WEP) 
and the strong (SEP). The WEP pertains to non-gravitational contributions to mass: namely. 
Standard Model contributions of nuclear and electromagnetic energy, plus quark masses and their 
kinetic energies. Nucleons of differing fractional electro-weak and nuclear binding energies might 
exhibit different couplings to gravity in the case of a WEP violation. The SEP extends the WEP 
to include gravitational self-energy of a. body, addressing the question of how gravity pulls on iiseif 
and. therefore, accessing the non-linear aspect of gravity. 

While the EP must hold true in GR, nearly all alternative theories of gravity predict a violation 
of the EP at some level. Efforts to formulate a quantum description of gravity genetically introduce 
new scalar or vector fields that violate the EP [14, 12] . These violations manifest themselves in the 
equations of motion for massive self-gravitating bodies, as well as preferred frame and preferred- 
location effects on the gravitational constant. GR may be the only metric theory of gravity that 
is dependent on the SEP holding true [68], distinguishing it from all other theories of gravity. 
Therefore, probing the validity of the EP is one of the strongest ways of testing GR. This test 
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is often considered one of the most powerful ways to search for new physics beyond the standard 
model [13]. 

Precision tests of the EP generally test the Universality of Free Fall (UFF): all test bodies 
have the same gravitational acceleration in a uniform gravitational field. Tests of the UFF are 
performed by comparing the gravitational accelerations ai and ct 2 of different test bodies: 


A a __ ai — <22 / Mq\ ( Mq \ 
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where Mq is the gravitational mass and A 1 1 is the inertial mass of each test body. Laboratory 
masses lack measurable gravitational self-energy, so the classical Eotvos type experiments which 
compare the acceleration of test bodies with different compositions only probe the WEI ’ [57], 

In the late 1960's. Nordtvedt recognized that the SEP could be tested by comparing the grav- 
itational acceleration of two massive bodies [48, 17], For each body, the gravitational to inertial 


mass ratio can be written as: 
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where U is the gravitational self-energy of the test body: 
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Me 2 is the body’s total mass energy, and rj is a dimensionless constant that is identically zero if 
the EP holds true. 

For a uniform sphere of radius R, U /M c 2 = — 3GM /5 Rc 2 , however due to their complex interior 
structure, the gravitational self-energy for astronomical bodies is generally computed numerically. 
An Earth model based on ;24l yields a self-energy of 171;: 
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while a Moon with a 20% iron core yields [72]: 
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Nordtvedt realized that a violation of the EP would cause the Earth and Moon to fall at different 
rates toward the Sun resulting in a polarization of the lunar orbit [44]. This polarization shows 
up in LLR as a displacement along the Earth-Sun line with a 29.53 day synodic period. Detailed 
solutions to the equations of motion for the Earth-Moon-Sun system [50. 15] find that the radial 
perturbation of the Earth-Moon distance due to an EP violation is 


Sr = -2.9427 x 10 10 
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where D is the angle between the mean longitude of the Moon and the mean longitude of the Sun 
as observed from the Earth (synodic period). Combining Eq. (2) with the estimated self-energy 
for the Earth and Moon, we find that Elq. (6) becomes: 

Sr = 13.1/; cos D [in], (7) 

Recent solutions using LLR data yield an HP test numerically comparable with present laboratory 
limits, at a part in 10 13 [74. 57]. Since the Earth and Moon not only have different, gravitational 



self-energy, but also have different compositions the LLR measurements alone do not provide a 
pure test, of SEP [49j . To separate the WEP and SEP effects and eliminate the possibility of 
a conspiratorial cancellation, the Eot-Wash group at the University of Washington performed a 
torsion balance experiment using test masses of similar composition to the Earth and Moon [4], 
Combining the torsion balance results with the latest LLR, analyses produced the best test, of the 
SEP to date (74, 75]: 

a(^) = (—2.0 ± 2.0) x 10~ 13 , (8) 

q = (4.4 ± 4.5) x 10 \ (Of 

Because Earth’s self-energy contributes 4.5 x 10~ 10 of its total mass, this translates to a SEP test 
of 0.04%. 


3 Variation of the strength of gravity 


In GR the coupling strength of gravity is taken as a constant: G. however, not all theories require 
such a constraint. Within scalar-tensor theories i he gravitational coupling can become a fund ion 
of a dynamical scalar field. The variation of the gravitational constant then depends on the 
cosmological evolution of the scalar field. The exact form of the variation will depend on the 
specific cosmological scenario being proposed [67, 26, 43]. 

A time variation of G will show up as an anomalous evolution of the orbital periods of astro- 
nomical bodies. This is easily seen from the Newtonian limit plus post Newtonian corrections to 
the two body orbital period: 
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where t is the angular momentum per unit mass, e is the orbital eccentricity, and m is a theory 
dependent, mass parameter of the two bodies [46], Taking the time derivative and ignoring the 
higher order terms we find 


For solar system bodies, we can safely ignore the mass term, except for a small rare of mass loss 
by the Sun (for compact objects like neutron st ars this term becomes important, [46j ). The ( terms 
account for the effects of torque acting on the orbit and must be carefully separated from the 
possible time evolution of G. 

The radial size of the Moon’s orbit is also affected by a varying G: 


7 


r 


G 

G' 


( 12 ) 


Both tidal-friction and a changing G influence the semimajor axis. However, one can separate 
the effects by taking into account their different proportional contributions to the orbital period. 
Combining Eqs. (11) and (12) one can solve for the angular momentum in terms of quantities 
measurable by LLR: 


which can then be used to calculate G/G. This simple analysis leads directly to an evolving range 
linear with time that was used for analysis of the initial years of LLR data. However, a changing 
G affects both the monthly lunar orbit and the annual Earth-Moon orbit around the Sun. Solar 



perturbations on the lunar orbit are also large. Secular change in the annual orbital period from a 
varying G accumulates as an orbital longitude perturbation evolving quadratically with time [51 . 
This t 2 effect on the phase of the solar perturbations provides a strong limit when measured over 
decades. 

Recent analysis of LLR data by the JPL group sets a limit on GjG = (f> ± 7) x 1. 0 — 1 ; year 

[65]. Similarly, Muller et al. find G/G = (2 ± 7) x 10 Ki /year and (7/(7 = (4 r 5) x 10 ' V/vear 2 

[35]. These limits translate to less than a. lf< variation of G over the 13.7 billion year age of the 
universe. Because the dominant effect for a variation in G is quadratic in time, continued LLR 
measurements will significantly improve this limit,. Additionally, a more optimal measurement 
schedule throughout the lunar month, now possible with APOLLO, will also put better constraints 
on a possible time variation of G [52]. 


4 Inverse square law 

The inverse square law (ISL) of gravity has been meaningfully tested over length scales spanning 
20 orders of magnitude, eliminating Yukawa-like couplings competitive with the strength of gravity 
from 10~ 4 to 10 16 meter length scales. The deepest probe of the ISL is from LLR. at a scale of 
~ 10® meters, where any new force must be weaker than gravity by more than ten orders-of- 
magnitude [37] . Short-range tests of the ISL have recently been prompted by the energy scale of 
the cosmological acceleration which suggest new-physics below 1 mm 111. 

Modern tests of Newton's inverse-square law of gravity often search tor an additional Yukawa 
contribution to the gravitational potential: 

V' (r) = -G ™ -St ( 1 4- ae .~'' /x | . (14 ' 

where a is the dimensionless strength and A is the. length scale. Such a potential would generate 
a precession of the Moon's perigee with frequency Sw j 1 j: 


where a is the mean radius of the Moon’s orbit. The agreement of geodetic prec ession with GR 
described below leads to a limit on an anomalous precession of 5io/u) < 1.6 x 10 _l1 . This translates 
into a limit on the strength of a new Yukawa potential of a < 5.9 x 10 -11 at A — a/2 where the 
lunar test is most sensitive. 

Recent analysis of LLR data includes specifically fitting for Yukawa perturbation terms in the 
equations of motion leading to a measurement of a = (3 ± 2) x 10" 11 at A = 4 x 10° km. While 
intriguing, this possible non-null result has yet to be thoroughly investigated 137 . 

5 The nature of spacetime 

The recent and unexpected measurement of the accelerating expansion of the universe has provided 
new motivation for exploring the nature of spacetime. Models that predict modification of gravity 
at large distances, such as bra.rt.e- world models, have recently become of interest 17}. These theories 
exhibit a strong coupling phenomenon that makes the gravitational force source dependent. These 
theories become testable at shorter distances where the coupling sets in for lighter sources ,lfy. 
The Earth-Moon system provides a testbed for investigating the nature of spacetime at solar- 
system scales. For example, GR predicts that, a gyroscope moving through curved spacetime 
will precess with respect to a rest, frame. This is referred to as geodetic or de Sitter precession. 
The Earth-Moon system behaves as a gyroscope with a predicted geodetic precession of 19.2 



msec/year. This is observed using LLR by measuring the lunar perigee precession [9]. The current 
limit on the deviation of the geodetic procession from the general relativity prediction is: K gp = 
(1.9 ±6.4) x 10~ 3 [741. This measurement, can also be used to set a limit on a possible cosmological 
constant: A < 10 -2b km~ 2 [58], which lias implications for our understanding of dark energy. 

It is also useful to look at violations of GR in the context of metric theories of gravity. Pa- 
rameterized Post-Newtonian (PPN) formalism provides a convenient way to describe a class of 
deviations from GR [53]. The most often considered PPN parameters are and 3: 7 indicates 
how much spacetime curvature is produced per unit mass, while 0 indicates how nonlinear gravity 
is (self-interaction). 7 and / 3 are identically one in GR. Also of interest are the preferred-frame 
parameters on and c*2 which are identically zero in GR [69. 54], 

Limits on 7 can be set from geodetic precession measurements [53], but the best limits presently 
come from measurements of the gravitational time delay of light, i.e. the Shapiro effect [59]. 
Doppler measurements to the Cassini spacecraft set the current limit on 7: (7 - 1) = (2.1 ± 2.3) x 
1CT 5 [8j. 

The Equivalence Principle parameter f) depends on the PPN parameters 3 and 7: 

/,^4/J — 7-3. (16) 


Combining the Cassini value for 7 with the LLR value for t) provides the best limit on 0: (0 — 1) = 
(I.2± 1.1 ) x 10 * [74]. Scalar tensor theories with attractors for the cosmic background scalar-field 
dynamics predict a residual 7 — 1 and perhaps 0 — 1 of order i0“'-10~ 5 today [12|, within reach 
of advanced LLR and spacecraft time-delay measurements. 

A non-zero preferred frame would show up as an oscillation of the lunar range at the sum and 
the difference of the anomalistic frequency and the annual period [10. 301. Recent analysis of LLR 
data sets the current limit on the PPN parameter ay ~ (—7 zt 9) x 10~ & [37]. LLR has also been 
used to set a limit 011 ao = (1.8 ± 2.5) x ICR' 3 [37], however, the close solar spin axis alignment 
with the total solar system angular momentum produces a much tighter constraint, on <*2 of order 
1CT 7 [45], 

Lunar laser ranging also places limits on the gravitomagnetic interaction, the same physical 
interaction that leads to the Lense-Thirring and Schiff precession phenomena as tested by precession 
of the LAGEOS orbital plane and by the precession 01 a gyroscope in Gravity Probe-B respectively 
[41 j . In the case of the lunar orbit, rotation is not involved, but rather translation of t he Earth and 
Moon point-masses in the solar system barycenter frame that produce 6 meter amplitude range 
signatures at both the synodic frequency and twice the synodic frequency. The amplitudes of these 
signatures are frame-dependent , reflecting the deep connection gravitomagnetism has with the 
covariant property of relativistic dynamics. Soffel et. al. showed the need for the gravitomagnetic 
term in the LLR equations of motion at the level of 0.15%. whether confined to a PPN context or 
allowed to vary independently [64]. If another experiment claimed a gravitomagnetic. or "frame- 
dragging" departure from GR at even the 1% level. LLR data would stand in conflict. [40j. 


6 Modeling and Analysis 

The analysis of LLR data requires a sophisticated model of the solar system eplienieris that also 
includes all the significant effects that contribute to the range between the Earth stations and the 
lunar retroreflectors. These models compute a range prediction and the partial derivatives of range 
with respect to each model parameter at the epoch of each normal point. The model predictions 
take into account orbital parameters, attraction to the Sun and planets, relativistic corrections, 
as well as tidal distortions, plate motion, and other phenomena that affect the position of the 
retroreflector and ground station relative to the centers of mass of the Earth and Moon [75] . Some 
of these parameters are measured by other means, but most are estimated from the LLR data. 



The range measurements are corrected for atmospheric delay and a weighted least square analysis 
is performed to estimate the ~ 170 parameters in the model [37, 34], most of which are initial 
conditions and masses of solar system bodies. LLR data is often combined with other spacecraft 
and planetary tracking data to further constrain the estimates or remove degeneracies. 

A number of models have been developed over the past 40 years. The model developed at 
the Jet Propulsion Laboratory (JPL) was one of the first programs for LLR analysis |70j and 
continues to be updated. It was recently used to produce limits on the SEP violation, time 
variation of the gravitational constant, and interior structure of the Moon 1 75, 7(i.. The open- 
source Planetary Ephemeris Program (PEP) is undergoing a major upgrade for LLR analysis at. 
the Harvard-Smithsonian Center for Astrophysics (CfA). It was also used for one of the first LLR 
analysis to test the SEP [60] and was recently used to test for Lorentz violation using LLR data 
[6]. A model developed at the University of Hannover was also recently used to produce limits on 
the relativity parameters, including the preferred frame PPN parameters [34, 37]. 

With the now routine operation of APOLLO at Apache Point, millimeter level data is being 
produced [5]. Unfortunately, none of the ephemeris models is currently capable of handling mil- 
limeter class data to maximum advantage [77]. New effects and error reduction techniques need to 
be added to the analysis tools that become important at the millimeter level. To take advantage of 
the next generation of LLR instruments, these codes will need to be further modified and rigorous 
theoretical work will need to be performed to permit tests of new ideas in physics. Substantial 
effort is also required to address the multitude of effects that will contribute at that level. 

Many of these effects will be scientifically interesting in their own right. In particular, the 
lunar interior models require significant improvement. There are also additional Earth effects that 
need further model development, such as the loading of the lithosphere by the atmosphere arid 
ocean which causes the observing station to move vertically (and horizontally) with the tides and 
weather. Models of these effects are available that are deemed accurate to better than 0.1 mm 
and tested in VLBI analysis software at CfA 61]. but need to be incorporated into the analysis 
programs. As ranging precision is further improved, more sophisticated atmospheric models and 
auxiliary measurements will need to be developed. Other important effects for advanced LLR 
analysis include solar radiation pressure, thermal cycling of the reflectors, solar tides on the moon, 
and solar mass changes [77;. 

7 Next Generation Lunar Laser Ranging 

The current lunar retroreflectors all lie within 26 degrees latitude of the lunar equator, and the 
most useful ones within 24 degrees longitude of the sub-earth meridian as shown in Figure 4. 
This clustering is sub-optimal, particularly with respect to measuring the lunar librations. In 
addition, the active LLR ground stations do not cover a large range of latitudes, further weakening 
the geometric strength of the observations. Additional observatories could improve the situation 
somewhat, but Mt. Strom lo in Australia is the only station capable of ranging to the Moon not 
situated at similar northern latitudes. Unfortunately, Mt. Stromlo has not been active in ranging 
to the Moon. 

The frequency and quality of observations also varies greatly with the facility and power of the 
laser employed. Most ranging over the recent past has occurred between three ground stations 
(MLRS, OCA, and Apache Point,) and one reflector (Apollo 15). The solar noise background and 
other issues make ranging to some reflectors possible only around the quarter-moon phase for most 
stations other than Apache Point, which has very good distribution among the reflectors. APOLLO 
is capable of ranging during all lunar phases, but it must share time on its 3.5 meter telescope 
with other programs. 

Improvements in the geometric coverage, both on Earth and on the Moon, will have a direct 
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impact on the science gained through LLR. Studies of the structure and composition of the interior 
require measurements of the lunar librations, while tests of GR require the position of the lunar- 
center of mass. In all, six degrees of freedom are required to constrain the geometry of the Earth- 
Moon system (in addition to Earth orientation). A single ranging station and reflector is insufficient 
to accurately determine all six degrees of freedom, even given the rotation of the Earth with respect 
to the Moon. The addition of one or more reflectors and one or more additional ranging stations 
in the Earth's southern hemisphere would strengthen the geometric coverage arid increase the 
sensitivity to lunar motion by as much as a factor of 4 at the same level of ranging precision [33] . 



Figure 4: Location of the lunar retroreflector arrays. The three Apollo arrays are labeled AP and 
the two Luna arrays are labeled LUN. ORI and SHK show' the potential location of two additional 
sites that would aid in strengthening the geometric coverage. 

Satellite Laser Ranging (SLR) began in 1964 at NASA’s Goddard Space Flight Center. Since 
then it has grown into a global effort, represented by the International Laser Ranging Service 
(ILRS) [55], The ILRS includes ranging to Earth orbiting artificial satellites, ranging to the lunar 
reflectors, and is actively working toward supporting asynchronous planetary transponder ranging. 
The current SLR network consists of over 40 stations worldwide, funded and operated by research 
organizations, universities, ' and foreign governments. It varies widely in its designs for tracking- 
mounts, telescopes, receivers, timing electronics, and laser transmitters. 

The development of NASA's Next Generation Satellite Laser Ranging (NGSLR) system has 
moved system operation from the regime of high laser energy/low repetition rate to high repetition 
rate single photon detection [21, 32J. This technique uses receive photons more efficiently and 
because of the higher return rate, minimizes acquisition time and enables closed loop tracking. 
The current laser in use (300 picosecond pulse width) limits single shot measurement accuracy to 
2 to 3 cm but because of the high return rate normal point data can be reduced to the millimeter 
level. As many as 12 NGSLR. stations are expected to be built and deployed around the world 
under NASA’s Earth Science Program in the coming decade. 

The development of this new network of SLR stations provides a potential opportunity to ex- 
pand the number of LLR stations. To be lunar capable, the SLR stations would either need to be 
upgraded with higher power lasers or new high cross-section retroreflectors and/or laser transpon- 
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Figure 5: Worldwide distribution of satellite laser ranging stations, participating in the Interna- 
tional Laser Ranging Service [55]. 

ders would need to be put on the Moon. "Faking advantage of the existing SLR infrastructure is 
a very compelling way to increase both the spatial and temporal LLR coverage at minimal cost., 
and would ensure continuous availability of LLR data through the indefinite future as it would not 
rely on unique facilities and individual investigators continuing operations. 

8 Next Generation Lunar Retroreflectors 

Remarkably, four of the five lunar retroreflector arrays are still visible and producing useful data 
after 40 years of exposure to the lunar environment. During that time, the precision of the range 
measurements was improved each time the ground stations were upgraded to the most advanced 
ranging technology at the time. Incredibly, it has taken nearly 40 years for the ground stations to 
catch up with the potential capability of the retroreflector arrays. 

The first; LLR measurements had a precision of about 20 cm. Since 1969. several stations have 
successfully ranged to the lunar retroreflectors and have increased the range accuracy by a factor 
of 10 to the level of a, few centimeters. Poor detection rates have historically limited LLR precision 
(not every laser pulse sent to the Moon results in a detected return photon). However, the rela- 
tively new APOLLO system uses the large collecting area of the Apache Point telescope, a tightly 
collimated beam, good atmospheric image quality, and has very efficient avalanche photodiode 
arrays such that thousands of detections are recorded (even multiple detections per pulse) leading 
to a statistical uncertainty of about, i mm for timescales of less than 10 minutes [5. 39 j . 

The dominant random uncertainty per photon received by APOLLO sterns from the physical 


11 



size of the arrays and their changing orientation due to the lunar librations. The incoming pulse 
from APOLLO will illuminate an entire array, but only one (sometimes a few) of the photons will 
be detected upon return. APOLLO cannot determine what area of the array contributed most of 
the returned light, so the tilt of the array with respect to the Earth spreads out the distribution of 
laser pulses. The typical array dimension of 0.5 m and a typical li brat ion angle of 6° translates to 
a full- width pulse-spread of about 0.1 m, or 670 ps in the round-trip time. As the Moon librates. 
the amount of spreading changes since the array is also changing its orientation with respect to 
the ground station. 



time offset (ns) 


Figure 6: Example data from APOLLO from the Apollo 15 array on November 19. 2007, in which 
6624 photons were collected in a 5000 shot run. The raw time events are shown in the top plot (with 
the initial predicted round trip time subtracted). The bottom plot shows the distribution of the 
outgoing pulses, which when convoluted with the relrorefiector tilt is consistent with tin measured 
returns shown in the central plot. The trapezoidal overplot represents the temporal spread due to 
the oriental ion of t he ret rorefleetor at t he time of t he observe t ion 39 j. Figure courtesy o: 1 homas 
Murphy. 

Modest improvements in the ranging technology will not significantly improve the range pre- 
cision as the array tilt will continue to dominate the error budget for the foreseeable future. In 
addition, new arrays with more (or less) cubes of the same size would result in no gain: doubling 
the physical dimension doubles the random uncertainty requiring four times as many photons, ex- 
actly what doubling the linear array dimension provides. Likewise, the reduction in return photons 
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would eliminate any benefit of going with a smaller array. 

To maintain the advantage that multiple cubes provide in response, but eliminate the issues 
with orientation, one can separate the cubes far enough that their responses do not overlap when 
seen by the Earth stations. A range separation of about 10 cm between the cubes should be 
sufficient to distinguish them in the data. Since the typical (for APOLLO) 5 microradian laser 
beam covers a 2 km spot on the moon, any reasonable spacing will result in illumination of the 
entire set of reflectors at once. The cubes could be coarsely surveyed individually to provide enough 
information to be acquired by the lower energy Earth laser stations, or the initial acquisition from 
the ground could be performed with the higher laser energy / larger telescope lunar laser ranging 
Earth stations which will have good signal-to-noise ratios. 

Large single cube corners can also be made to provide similar return rates as the Apollo arrays 
without significant pulse spreading. The response from a 7.6 cm cube would be 16 times larger than 
that of the Apollo 3.8 cm cubes. However, simply making solid cubes larger increases their weight 
by the ratio of the diameter cubed. The additional size also adds to thermal noise effects and 
decreases the cube’s divergence: a very narrow divergence can cause the return spot to completely 
miss the station due to velocity aberration. Spoiling (making the dihedral angles of the cube 
different from 90 degrees) can compensate for the velocity aberration but reduces the effective 
cross section. 

Solid cube corner retroreflectors (up to 11 cm) have flown on over a hundred missions, for both 
satellites and lunar laser ranging. Recent tests of the 10 cm cube shown in Fig. 7 has demonstrated 
it meets relevant requirements for the lunar environment [11]. Designs for the housing are still in 
development [22], 



Figure 7: A 10 cm solid cube corner reflector was recently qualified for the lunai environineni 
Also shown for comparison is a 3.8 cm Apollo engineering model cube corner 1 1 . Photo courtesy 
of Douglas Currie. 

The main disadvantage of large solid cubes is that the thermal requirements become very 
challenging because of the temperature dependence of the cube material’s index of refraction. A 
promising alternative is to use hollow (open) cube corners. Since hollow cubes are reflective, the 



index of refraction problem goes away. They also potentially weigh less, have smaller t hernial 
distortions, and do not introduce significant polarization effects. They, therefore, can be made 
larger without sacrificing as much in optical performance. Hollow cubes have flown on a few 
missions, but are generally not used on satellites for laser ranging because of a lack of test data 
and some indications of instabilities at high temperatures. A recent program at NASA Goddard 
Space Flight Center is looking at applying advanced bonding techniques for space optics that have 
the potential for mitigating these problems. 

Isolation from ground motion and thermal changes are also key for going beyond the Apollo 
array capabilities. Each reflector should be rigidly grounded to directly sense lunar body motion 
and be located far enough away from normal human activity to avoid vibration and contaminat ion 
(dust) from affecting the cubes. To provide thermal stability, the retroreflectors could be thermally 
coupled to the ground below the surface layer. As shown in Fig. 8, measurements from the Apollo 
15 and 17 heat flow probes indicate that the large diurnal temperature fluctuations are negligibly 
small at depths below 0.8 meters [31, 27], To take advantage of this stability, one could drill a hole 
about a meter deep and insert a rod with high conductivity, low coefficient of thermal expansion 
to anchor the retroreflector which would be mounted to the exposed end of the rod. A thermal 
blanket positioned over the lunar surface around and below the retroreflector would also reduce 
the thermal fluctuations induced from the surrounding regolith. 



Figure 8: Temperature fluctuations in the lunar regolith as a function of depth from Apollo 15 and 
17 measurements. Hatched areas show day-night temperature fluctuations. Below about 50 cm 
there was no observable temperature fluctuation due to the lunar day- night temperature cycles. 
[31]. 


9 Laser Transponders 

Laser transponders are currently being developed for satellite laser ranging, but they can also be 
deployed on the lunar surface as an alternative to retroreflectors [20 j. Transponders are devices 
that both send and receive predictable signals and can be used for ranging and time transfer. They 
have approximately a r 2 link advantage over direct ranging loss of 1/r 4 , as the signal is propagat- 
ing in only one direction before being regenerated. With the development and inclusion of laser 
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communications for spaceflight missions, it is logical to include an optical transponder that uses 
the same opto-mechanical infrastructure with minimal impact on the mission resources. These 
instruments could be used to support the relativistic and lunar science in addition to providing 
communications support to astronauts and/or other scientific instruments. These lunar instru- 
ments would also provide a pathfinder for applications on Mars and other planetary bodies where 
the use of passive retroreflectors is not practicle. 

Asynchronous Laser Transponders (ALT) are relatively simple devises that have good potential 
for space applications. The ground and remote stations of an ALT fire their lasers independently 
(as opposed to an echo transponder which works by sending back a timing signal with a fixed 
delay from the receipt of the base-station signal). This allows the use of free-running lasers on 
the space-based terminal to operate at their most efficient repetition rates. It is also simpler and 
potentially more reliable than other types of laser transponders. It does, however, require a clock 
on both ends. 



Figure 9: An Asynchronous Laser Transponder is under development at NASA Goddard Space 
Flight Center that is compatible with the Next Generation Satellite Laser Ranging system 1 9 


Several interplanetary laser transponder experiments were successfully performed from the 
NASA Goddard Space Flight Center satellite laser ranging facility. The first (two-way. 24 x 10° km 
link) utilized the Mercury Laser Altimeter on i he Messenger spacecraft |63j and the second (one- 
way. 80 x 10® km link) utilized the Mars Orbiting Laser Altimeter on the Mars Orbiter spacecraft 
1781. More recently, two-way ranging was successfully performed using the Lunar Obiter Laser 
Altimeter' (LOLA) on the Lunar Reconnaissance Orbiter (LRO) orbiting the Moon [32;. These 
experiments have proven the concept of being able to point both transceivers, detect the photons, 
and retrieve useful parameters at low-link margins over interplanetary distances. 

An ALT, shown in Fig. 9, is currently under development at NASA Goddard Space Flight 
Center that uses technology derived from the Next Generation Satellite Laser Ranging (NGSLR) 
system [19], Other efforts at Goddard and other institutions are currently underway that combine 
laser ranging and laser communications. The robust link budget combined with the potential 
compatibility with NGSLR would open up to the possibility of using a large number of existing 
ground stations, which would not only increase the scientific potential, but also significantly reduce 
the ground station and operational costs. 


10 Outlook 

In addition to the tests of gravity, new LLR instruments would be of great benefit to studies of t lie 
interior structure of the Moon. Tighter constraints resulting from more complete tracking could 


Figure 10: NASA’s Next Generation Satellite Laser Ranging (NGSLR) system is being used to 
laser range to the Lunar Reconnaissance Orbiter (LR.O) orbiting the moon [32]. 

aid in the search for a solid inner core. The second-degree tidal lunar Love numbers are detected 
by LLR. as well as their phase shifts. From past measurements, a fluid core of ~ 20% of i he 
Moon’s radius is indicated. A lunar tidal dissipat ion of Q - 30 zr 4 has been reported to have a 
weak dependence on tidal frequency [73;. Evidence for the oblateness of the lunar litiid-core/solid- 
mantle boundary may be reflected in a century-scale polar wobble frequency. The lunar vertical 
and horizontal elastic tidal displacement Love numbers hi and li are known to no better than 
25% of their values, and the lunar dissipation factor Q and the gravitational potential tidal Love 
number no better than 11%. These values have been inverted jointly for structure and density 
of the core [29, 28], implying a liquid core and regions of partial melt in the lunar mantle. 

Lunar interior studies have arguably suffered the most from the clustering of the Apollo arrays 
on the central portion of the moon. Placing retroreflectors far from the Apollo arrays, at a pole 
or a limb, would improve the measurements by up to a factor of 4 at the same level of ranging 
precision as is currently performed [33]. 

The advancement of active laser ranging systems also opens up the possibility of precision 
ranging beyond the Moon. Laser ranging to Mars can be used to measure the gravitational time 
delay as Mars passes behind the Sun relative to the Earth. With 1 cm precision ranging, the PPN 
parameter 7 can be measured to about 10~ 6 . ten times better than the Cassini result [66 i. The 
Strong Equivalence Principle polarization effect is about 100 times larger for Earth-Mars orbits 
than for t he lunar orbit. With 1 cm precision ranging, the Nordtvedt parameter. // = }/? - 7 - 3. 
can be measured to between 6 x 10 -b and 2 x 10~° for observations ranging between one and 
ten years [2]. Combined with the time delay measurements this leads to a measurement of PPN 
parameter 3 to the 10 -6 level. Mars ranging can also be used in combination with lunar ranging 
to get more accurate limits on the time variation of the gravitational constant. 

The ephemeris of Mars itself is known to meters in plane, but hundreds of meters out-of-plane 
[30] . Laser ranging would get an order of magnitude better estimate, significant for interplanetary 
navigation. Better measurements of Mars’ rotational dynamics could provide estimates of the size 
of its core [25], Mars’ elastic tidal Love number is predicted to be less than 10 cm, within reach 
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of laser ranging. There is also an unexplained low value of Q. inferred from the secular decay of 
Phobos’ orbit that is a constraint, to the present thermal state of the Mars interior 10 . Laser 
ranging to Phobos would help solve this mystery. 

LLR remains one of the best tests of gravity in the weak field and promises to cont inue to be 
a key tool for many years to come. Four of the five lunar retrorefiectors remain visible today and 
continue to produce valuable data. Advances in ranging technology has finally reached the point 
where the precision of the data is being limited by the physical characteristics of the lunar arrays. 
The natural next step in LLR is to place new retrorefiectors and/or laser transponders on the Moon 
at sites far from the Apollo arrays that have a high enough return rate to take advantage of the 
SLR network of ground stations. With the retrorefiectors and transponder technology available 
today, these new instruments could easily support laser ranging and advances in ground station 
technology for another productive 40 years of LLR. 
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